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If  thermal  conductivity  and  specific  heat  are  taken  as  linear  functions  of 
temperature,  a  nonlinear  heat-conduction  equation  results.  For  small  non- 
linearities  an  approximate  first-order  analytical  solution  may  be  obtained 
in  certain  ccses.  The  present  analysis  deals  with  one-dimensional  prob¬ 
lems  with  periodic  boundary  conditions.  Only  the  steady-state  solution 
(i.e,,  one  which  is  periodic  in  time)  is  considered.  Solutions  are  obtained 
for  the  following  cases:  1)  Semi-infinite  solid  with  sinusoidal  boundary 
temperature,  2)  thick  slab  with  sinusoidal  temperature  at  one  boundary 
and  constant  temperature  at  the  other,  and  3)  thick  slab  with  prescribed 
heat  flux  (a  constant  term  plus  a  sinusoidal  term)  at  one  boundary,  con¬ 
stant  temperature  a*  the  other.  The  effects  of  the  nonlinearities  are  dis¬ 
cussed;  they  are  found  to  be  surprisingly  small. 
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*  1 

One-Dimensional  Quasilinear  Heat  Flow 
With  Boundaiy  Conditions  Periodic  in  Time  >m  $.  Thompson 

NOMENCLATURE 

The  following  nomenclature  is  used  in  this  paper: 

a  =  Temperature  coefficient  for  k  in  Equation  [s] 

=  Amplittides  of  constant  and  first  harmonic  terms  of 
*  heat  flux  input.  Equation  £43-3 

b  =  Temperature  coefficient  for  j9  in  Equation  16] 

c  =  Specific  heat  (for  unit  mass) 

f  =  Time-independent  term  in  temperature  solution 

k  -  Thermal  conductivity 

*2t0  =  Value  of  k  at  reference  temperature  (assumed  zero) 

=  Thickness  of  slab 
=  Heat  flux 
^  =  Amplitude  of  n—  harmonic  of  temperature 
m  =  eg  =  .Volumetric  specific  he«t 

=  Value  of  £  at  reference  temperature  (assumed  zero) 

Tq  =  Amplitude  of  temperature  input.  Equation  [171 

tV-  =  /  VT  ICq  y,  ,  =  Amplitude  of  linear  solution  at  boundary 

T^  =  Constant  temperature  at  x  =  L* 

t  =  Temperature 

u  -  Transformed  temperature,  defined  by  Equations  W  and  [7] 
x  =  Distance  coordinate 

=  k  /  eg  =  Thermal  diffusivity 
in  =  nm/2  «„ 
i  =  (b-a)  Bq  /  2k„  y, 

|°  =  Density 

IT  =  Time 
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4>  =  Coefficient  of  cosine  term  of  n—  temperature  harmoniA  NT,S  cra&i 

r  n  ;  I  DTIC  tab 

Unannounced 


th 


*X  n  =  Phase  angle  of  nr1*  harmonic  in  temperature  solution 
Y  x\  ~  Coefficient  of  sine  term  of  n —  temperature  harmonic 
***  =  Angular  frequency 

INTRODUCTION 

The  general  equation  for  heat  conduction  in  an 
isotropic  material,  without  heat  sources  or  sinks,  is 
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div  (k  grad  t 


where  t  is  temperature,  \  is  time,  k  is  thermal 
conductivity,  and  si  =  c  p  the  volumetric  specific  heat, 
i.e,,  the  product  of  specific  heat  and  density.  In  the  case 
where  k  and  s  are  constants  this  becomes 


'at 


where  o(S  k^  Is  a  constant.  This  is  a  linear  oartial 
differential  equation;  solutions  have  been  obtained  by 
standard  methods  for  many  cases  of  practical  interest  (see 
for  example  ^(1)).  However,  in  the  more  general  case,  both 
k  and  £  in  Equation  £l]  are  functions  of  temperature;  hence 
the  differential  equation  is  nonlinear. 

In  the  nonlinear  case  k  and  £  may  readily  be 
combined  into  a  single  temperature-dependent  quantity  through 
the  transformation  used  by  van  Dusen  (2) 
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\  } 

where  k  is  the  value  of  k  at  temperature  t,,.  It  follows 
that  k  grad  t  =  k^  grad  u  and  (3t/3t)»  (K/k)  (3u/33t) 
and  hence 

-  <“> vV  “  ft 

where  <*(u)  =  k  (u)/s  (u).  For  the  t ime- independent  case, 

\7  2  u  =  0  ;  thus,  after  the  boundary  conditions  have  been 
expressed  in  terms  of  u  ,  the  problem  reduces  to  the  familiar 
linear  one, 

In  the  transient  case  the  scope  of  available  analytical 
solutions  is  much  more  limited,  Hopkins  (3)  has  suggested  the 
method  of  successive  approximations  and  applied  it  to  a 
number  of  examples.  Another  analytical  technique  has  been 
the  Boltzmann  transformation  (U,  which  converts  the  nartial 
differential  equation  into  an  ordinary  one;  this  is 
limited  to  the  one-dimensional  case  with  rather  specialized 
Initial  condition  and  boundary  conditions.  Additional 
solutions  for  various  cases  are  given  by  Crank  (7)  and  by 
Friedmann  (5*6)  who  includes  an  extensive  bibliography.  Much 
of  the  work  in  the  field  is  based  on  numerical  solutions;  most 
of  it  limited  to  the  one-dimensional  problem. 

The  case  in  which  boundary  conditions  are  periodic  in 
time  has  not  received  as  much  attention.  However,  if  the 
initial  transient  is  neglected,  the  time-dependent  steady  state 
solution  can  readily  be  found,  under  physically  reasonable 
restrictions,  Vernotte  (8,9)  has  treated  a  semi-infinite 
bar  with  sinusoidal  heat  flux  input  and  has  briefly  discussed 
the  case  of  square-wave  heat  flux. 

The  present  naper,  based  on  the  fuller  analysis  by 
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one  of  the  authors  (10),  is  an  extension  and  amplification 
of  Vernotte’s  work*  It  treats  additional  problems,  presents 
more  explicit  results,  and  includes  second  harmonic  terms* 
(However,  it  deals  only  with  a  one-dimensional  solid  rather 
than  a  bar  with  convection  losses  at  the  cylindrical  surface.) 
Qualitative  conclusions  are  discussed  from  two  different 
view-points:  a)  the  case  in  which  the  thermal  properties 
are  known  and  the  temperature  solution  is  desired  (e.g*  the 
periodic  solar  heating  of  the  crust  of  the  earth  or  the 
periodic  heating  of  a  cylinder  wall  in  an  internal  combustion 
engine);  b)  the  case  in  which  the  required  temperatures  are 
experimentally  measured  in  order  to  determine  the  thermal 
properties  * 


PROBLEMS  AND  ASSUMPTIONS 

Three  problems  will  not  be  considered: 

a)  A  one-dimensional,  semi-infinite  solid  with  the 
surface  temperature  sinusoidal  in  time. 

b)  A  one-dimensional,  thick  slab  with  sinusoidal 
temperature  variation  at  one  boundary  and  constant  temperature 
at  the  other.  (The  term  "thick  slab"  implies  a  negligibly 
small  periodic  temperature  variation  near  the  second  boundary. 

c)  A  one-dimensional  slab  with  the  right  boundary  at 
constant  temperature.  At  the  left  boundary  there  is  a 
prescribed  heat  flux  input,  consisting  of  a  constant,  tern 
plus  a  sinusoidal  term.  This  example  may  be  useful  in  the 

experimental  measurement  of  thermal  conductivity. 

Two  basic  assumptions  are  made  throughout  the 


analysis 


1)  It  is  supposed  that  the  periodic  input  has 
been  applied  long  enough  so  that  the  transient  has 


died  out.  In  other  words,  only  the  steady-state, 
time-dependent  solution  is  obtained;  it  is  assumed 
that  the  solution  is  periodic  with  the  same  period 
as  involved  in  the  boundary  conditions.  Hence,  no 
initial  conditions  appear  in  any  of  the  cases 
considered. 

2)  It  is  assumed  that  the  temperature- 
dependent  r ••‘onerties  are  linear  functions  and  that 
the  variations  are  small.  ore  specifically,  the 
thermal  conduct! vity  k  and  the  volumetric  specific 


taken  as 

-  kQ  (1+2  at) 

[si 

=s  (1  +  2  bt) 

[6l 

where  kc,  s c,  a  and  b  are  constants.  Since  the 
variation  is  assumed  small,  2at <<1  and  2bt«l . 

Hence  terms  independent  of  a  and  b  are  considered 
as  nero-order  terms,  linear  terms  in  these  parameters 
are  considered  as  first-order  in  smallness,  quadratic 
terms  in  these  parameters  are  taken  as  second  order  and 

assigned  negligibly  small. 

METHOD  OP  ANALYSIS 

The  method  of  solution  is  that  of  "harmonic  balance." 
This  technique  is  the  same  as  that  employed  by  Vernotte  (8,9) 
except  that  he  deals  with  amplitude  and  phase  angle  for  each 


frequency  while  the  present  analysis  deals  with  amplitudes 
of  sine  and  cosine  terms,  "Harmonic  balance"  is  a  familiar 
technique  for  ordinary  nonlinear  differential  equations 
(e.g.,  references  (11)  through  (llj.). 

It  will  be  convenient  to  treat  the  problem  in  terms 
of  the  transformed  variable  u  defined  by  Equation  [3]- 
Using  Equation  [S]  for  k  in  this  transformation  yields 

u  =  t  +  at2  =  t  (1  +  at)  [7] 


If  second-order  terms  are  neglected  (as  stated  in  Assumption  2), 
1+au  *  1+at  and  the  Inverse  transformation  is 


t 


u  (1  +  au) 


-1 


u  (1  -  au) 


In  the  same  way.  Equations  [s]  and  [6]  become  k  "Sr  (1+2  au) 
and  s  =t  so  (1  +  2  bu).  Substituting  these  expressions  in 
the  one-dimensional  case  of  Equation  gives 


or 


[9] 


f*  ~  H  +  *  (*>-*)*  'll*  ,  J.  li*  9U* 

<*.  o(# 


[10] 


where  oC  s  k  /  s  • 

0—0—0 

As  stated  in  Assumption  1,  the  solution  will  be 
assumed  periodic  and  may  thus  be  written  as  a  Fourier 
series,  i.e. 

ia  (*  */  s  f  (*)  1>%  W  Cpi  W  -r  )►,  (X.;  w-r 

+  ^(%;*os  2wt  Sin  aw'C  O'] 
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CUT  higher  terms  are  included  in  the  series,  they  will  be 
found  to  be  at  least  second-order  for  the  problems  under 
consideration  here,)  Inserting  Equation  M  in 
Equation  [10]  now  yields  the  approximate  expression 

si*  2wt 


if  .  ^  c>  mlC  +  *,'n  wt  1*  C°S  2W-tr  T-T 

**  ^ 


JL  p  sin  wt  +  CAiUt  -  2fju  sin  2wt  -+2)^w  cos  2w^] 


£-  2"f  £  c*l  *n  w“t  +  2f  yt  u»  Cos  ttf't  -  if\-  f* )  2^^  +  2  ^  jf»(  cos  2u»  "tj 


(Since  the  last  bracket  is  multiplied  by  (b  -  a)  and  since 
and  fa  may  be  shown  to  be  first-order  terms,  their 
product  has  been  omitted  in  the  bracket,) 


Equating  coefficients  of  the  time-independent 
terms  and  of  the  respective  sine  and  cosine  terns  now  gives 
the  five  ordinary  differential  equations 


[12] 


g.f  (b-a,)  uj 

<*• 


% 


^ 

-  M 


2 

«<• 


Eis] 
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Solution  of  a  specific  problem  now  consists  in 


expressing  the  boundary  conditions  in  terms  cf  u,  comparing 
these  expressions  with  Equation  Tn]  to  determine  the  boundary 
conditions  on  f,  t, ,  1>i  ,  and  Yz  »  solving  the 

five  ordinary  differential  equations  (Equations  fisT  through 
[16]  )  subject  to  these  conditions  and,  finally,  expressing 

the  solution  in  terms  of  t  .  The  various  steps  will  be 
explained  in  some  detail  in  the  first  problem  below,  while 
the  other  cases  will  be  outlined  more  briefly* 


SOLUTIONS 

a)  Semi-Infinite  Solid,  Periodic  Surface  Temperature 

This  problem  is  mathematically  defined  by:  1)  the 
one -dimensional  form  of  Equation  [1]  ,  2)  the  expressions 
for  k  and  js  ,  Equations  [5]  and  [6]  ,  and  3)  the 
boundary  conditions 


t  [O'X)  *  T0  sin  W  *  |"17] 

tK*;  l  [l8] 

The  first  two  items  are  combined  in  Equation  [io]  ,  which 
in  turn  led  to  Equations  [l2]  through  [l6].  The  third  item, 
the  boundary  conditions,  may  be  transformed  by  Equation  p7 ^ 
to  become 


U  (o 


>  x)  =  4-  T0  Sm  W-t  _ 

/  •  a  2 


I  *  C00/  V I 
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By  comparing  u  (o,  X  )  with  Equation  [ll}  ,  it  follows  that 

fw- 


[19] 


*kW=T. 


[a] 


fz  (®i  *  o 


C21] 


u  (  )  indicates  that  f, 

remain  finite  at  infinity* 

The  solution  of  Equation 

'T.1 


f  (*■) 


/  r> ,  *k  “nd  n  a11 

[12]  is  obviously 


Since  f  (x)  itself  is  first-order  in  this  case,  the  right 
hand  sides  of  Equations  [13]  and  are  clearly  negligible. 

Hence  these  two.  equations  with  the  boundary  condition  given 
by  Equation  [20l  reduce  to  the  well-known  linear  case;  the 


solutions  are 

-r  * 

</>t  (*)  *  -  \  e 

[23] 

f,  *  T0€  *'  Coi 

*,* 

where 

[25] 

Inserting  Equations  [23]  and  [2l|]  in  [if>3  and 
£l63  yields 
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A 

.U* 

2W 

r» 

(b-a)(o  T0* 

-2*,* 
e  Sin 

s*,* 

A  + 

geo 

** 

__  (b-*)  to 

e  c%> 

2 

<»%* 

*• 

with  boundary  conditions  given  by  Equation  [21]  .  The 
solutions  are 


»4*  -  <1^3*  e'ilr,xc,  2,(IC 

%<*>  *  -  e'^'"  %*  -  tdi*  „*  2^ 


where 

X2  a  W^.  *  >15  *1 


The  solution  in  terms  of  u  is  now  obtained  by 
inserting  Equations  [ 22 ]  ,  [23]  ,  [2i+]  ,  [26]  ,  and  [27I 
in  Equation  [ll]  and  finally  converted  back  to  t  by  the 
inverse  transformation.  Equation  [a].  When  second -order 
terms  are  neglected,  the  result  is 


[26] 

[-27] 

[28] 


t2  -a*,*  -  *( 

2lLs  «  +  T0  e  s«n  Xk 


+  SllhLI0  Fe  '  c  o  i  (a  w  -t  -  a  y(  x)  -  <*  (s  x)l 

a  L  * J  £9] 

It  is  of  interest  to  consider  the  amplitude  Ibj 
and  the  phase  angle  of  the  second  harmonic  term*  In 

reference  (10)  it  is  shown  that  these  quantities  are  closely 
approximated  by 
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1.1 


Fig.  1  Temperature  solution  for  different  positions 
in  a  semi-infinite  solid  with  sinusoidal  boundary 
temperature,  given  by  Equation  (29),  for  the  case 
2(a)  T0  =  0.25,  b  =  0.  The  first  harmonic  is  iden 
tical  to  the  linear  solution. 


\ 

- R./-U 

\ 

- F./t. 

- R,/T. 

___  . 

*,X — - 


Fig.  2  Amplitudes  of  time -independent 
term  F0,  first  harmonic  Rj,  and  second  R2 
in  a  semi-infinite  solid  with  sinusoidal  boun¬ 
dary  temperature  for  die  case  2(a)  T0  =  0. 25. 
(b)  =  0. 


Rs=  AZ  (2«-t>)X2  Jr,*:  e'1'71 

-X  =  -  f  -  1-71  JT,X  +.03 


Prom  equation  on  it.  readily  follows  that  the  maximum 
valxie  of  R  occurs. at  =  .59  and  is  given 'by 


-2 

(R2^max 


.09  (2a-b)T 


Figures  1  to  3  shov/  various  aspects  of  the  nonlinear 
solution  for  the  case  in  which  2aT0  =  .25,  i.e.,  the 

thermal  conductivity  varies  by  +  ZScr-  about  its  mean  value. 
The  effects  of  the  nonlinearity  appear  considerably  smaller 
than  might  be  expected  from  a  preliminary  estimate. 


bj  Thick  Slab,  Periodic  Surface  Temperature 

The  problem  of  the  thick  slab  with  periodic  surface 
temperature  is  identical  with  the  previous  case  except  that 
Equation  [18]  is  replaced  by  a  boundary  condition  at  the 


right  hand  surface  X  -  L 

#  i*e*a 

t  (<-,*) -Tu 

[33] 

where  it  will  be  assumed 

1 ZL  I  <  I  !o  I  •  By  Equation  [7]  , 

this  becomes 

u  (L,  )  =  TL  +  aTL2 

[34] 

Comparison  with  Equation  [n]  now 

2 

f  (L)  =  TL  +  a  Tl 


Fig.  3  Maximum  and  minimum  temperatures  in 
a  semi-infinite  solid  with  periodic  boundary  tem¬ 
perature  for  die  case  2(a)  T0  =  0. 25,  (b)  =  0.  The 
linear  case  is  shown  for  comparison. 

Shows 

[35] 


($0-)  =  fcW-  (L)  (L)  =  O 


The  problem  now  reduces  to  solving  Equations  |l23  through 
[16)  with  boundary  conditions  given  by  Equations  [19],  [20], 
[21]  .  [35]  ,  and  [36] . 

The  time- independent  part  f  (x)  j.s  easily  found  to 
be 


since  the  term  "thick  slab"  will  be  taken  to  imply 


-8.X  -  Vax 

Hence  terms  in  £  or  e  may  be  assumed  zero  at 

the  right  boundary.  Furthermore  terms  in  (  Jf,  L)  1  may  be 
taken  as  first-order  in  smallness* 
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Inserting  Equation  [37]  in  Equation  [13]  now  yields 

4  -  ?  K  -  MW 


to] 


w4ll  be  approximately  equal  to  the  solution  of  the  linear 
problem  as  given  by  Equation  [2lJ .  Thus  the  order  of  magnitude 
of  the  last  factor  in  Equation  [39]  is 

The  maximum  value  of  the  term  on  the  right  is  .37  T  (  6,  L)”1; 

— o  ' 

hence  it  is  a  first-order  term.  Since  the  right  side  of 
Equation  [39]  also  contains  a  factor  of  (b-a)T^,  the  product 
is  negligible.  Similarly  the  right  hand  side  of  Equation  [ikl 


vanishes. 


Thus  the  remainder  of  the  problem  now  becomes 


identical  to  the  previous  case  and  the  final  temperature  solution 


t  (*  -t)  -  \  (1 + »\  -  *-^r)  (r) 


c  1  (w't-jrx) 


x_ar  .ay,*  -JL% 

v  J*  1  e  co*  e  cos  (awt-jj*)] 


Ro] 


c)  Thick  Slab,  Periodic  Heat  Flux. 

The  boundary  conditions  for  a  thick  slab  with  periodic 


heat  flux  will  be  taken  as 


!>,*>- -k(3y-  B-+B- 5,n 


Dal 
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[1*1 


(The  second  condition  is  chosen  to  make  u  vanish  at  the  origin,) 
Prom  Equation  [3  3  ,  it  follows  that  ]S  is*  (^/^)  j 

with  this  relation  and  Equation  [7]  ,  the  boundary  conditions 


become 


f2>*\  _  Ik  _  — 

(h.  ■  “  ■  ■ 


sin  to  "c 


'X*o  •  ' 

u  (L,  't  )  =  -  B  L/k 

o  o 

Hence  the  boundary  conditions  on  Equations  [12]  through  [16]  are 


Do] 

W 


(it  \  B. 

1  K  > 

f  (L).  -  %Jr 

[45] 

(^1  *  °  /  (^)o=  '  ^ 

JKW 

rw] 

ft?),*  0  “  (S).  / 

!*,  00  =  0  = 

147] 

The  solution  of  Equation  [i2^  is  now 

>T 

it 

1 

1  w 

mo 

0  1 

1  * 

[Ip] 

Substituting  this  result  in  Equation  £l3l  yields 

a1*  w 

1  «.  ed  '  *  L  w 

[49] 

ax' 


Since  the  coefficient  on  the  right  side  of  Equation  [1^91  is 
first-order  and  since  ft \  differs  from  the  linear  solution, 

t  by  first-order  quantities,  it  is  sufficiently  accurate 
to  use  Ufj  in  this  term  and  obtain 
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2(b-*)  B#u>  X 


When  Equations  and  [5X1  are  solved  with  boundary 

conditions  given  by  P*1  ,  the  result  is 

y;C*)  -  T.*e'‘r,Vi){D*sW+HW1]«- 

_  i  (»,*)*  *i«  (If.K--*-  ?•  ♦  a)} 


where 


T.»  =  B,  /V?  k. », 
ss  (b-»)  B. /i  k.», 

M 

is  the  amplitude  of  the  linear  solution  at  the  boundary 
while  $  is  a  dimensionless  parameter  giving  one  measure  of 
nonlinearity,  and  may  nowbe  found  as  before. 

In  terms  of  t  ,  the  final  result  becomes 


Fig.  4  Temperature  solution,  given  by  Equation 
(56),  at  g  ,  x  -  0. 5  in  a  thick  slab  with  input 
flux  consisting  of  a  constant  plus  a  sinusoidal  term 
for  the  case  2(a)  TQ  =  0.25,  (b)  =  0,  BQ  =  Bj.  The 
linear  solution  is  shown  for  comparison. 


term  F0,  first  harmonic  Rj,  and  second  har¬ 
monic  Rg  in  a  thick  slab  with  input  flux  con¬ 
sisting  of  a  constant  plus  a  sinusoidal  term 
for  the  case  2(a)  T0  =  0. 25,  (b)  =  0,  B0  = 


t  (x,Tr)  - 


o. 


T/ 

a 


+  X *  e~,,Xjp+§+  *(*,*)]  «■»  V?  %{*%)  sin 


+  T*  (2u,r.^+  [SB] 

9 

Some  aspects  of  this  solution  are  shown  in  Figures  4  and  5* 

CONCLUSIONS 

a)  The  method  presented  here  offers  a  rather  general 
approach  for  finding  the  steady-state  solution  when  boundary 
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conditions  are  periodic  and  variations  of  the  thermal  parameters 
are  small*  However,  the  algebra  may  become  quite  complicated, 
particularly  in  two-dimensional  and  three-dimensional  problems* 

b)  The  most  important  qualitative  effect  appears  to  be 
the  introduction  of  a  second  harmonic  term;  in  the  linear  case, 
such  a  term  does  not  appear  at  all.  Both  the  time-independent 
term  and  the  first  harmonic  may  also  be  modified;  in  fact,  in 
some  cases,  it  appears  that  a  time-independent  term  is 
introduced  by  nonlinearity.  However,  it  must  be  recalled  that 
the  zero  level  of  temperature  is  arbitrary;  hence  the 
time-independent  term  can  be  regarded  as  a  quantitative  change 
rather  than  a  qualitative  one, 

c)  The  quantitative  effects  of  the  nonlinearities  are 
surprisingly  small;  i^t  appears  that  the  thermal  properties  may 

safely  be  treated  as  constants  in  some  cases  where  their 
variation  is  quite  large.  For  example,  when  the  thermal 
conductivity  at  the  surface  varies  by  +  about  its  mean 
value,  the  maximum  amplitude  of  the  second  harmonic  may  be 
only  2 %  or  3/°  of  the  amplitude  of  temperature  variation  at  the 
surface*  (The  effect  on  the  time-independent  term  is  roughly 
twice  as  great;  in  the  last  problem  this  was  accentuated  by  the 
constant  term  in  the  heat  flux  input*)  The  explanation  seems 
to  be  that  the  conductivity  varies  more  or  less  symmetrically 
about  its  mean  value  kQ  ;  if  the  -variation  in  k  were  pro¬ 
portional  to  the  square  of  the  temperature,  nonlinear  effects 
might  be  considerably  greater* 

d)  In  general,  variation  in  thermal  conductivity  k 
has  more  effect  than  variation  in  the  volumetric  specific  heat 
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£  .  The  time-independent  part  of  the  solution,  if  any,  is 
unaffected  by  £  • 

e)  Since  the  second  harmonic  is  a  characteristically 
nonlinear  term  and  vanishes  for  the  linear  case,  experimental 
measurement  of  the  second  harmonic  amplitude  seems  to  provide  a 
sensitive  method  of  determining  the  coefficients  a  and  b  , 
which  describe  the  temperature  variation  of  k  and  £  , 
respectively.  However,  several  possible  drawbacks  should  be 
noted.  First,  it  may  be  difficult  to  generate  a  sinusoidal 
temperature  or  heat  flux  input  with  sufficient  accuracy. 

Secondly  on  the  basis  of  the  second  harmonic  alone  it  is 
d.fficult  or  impossible  (depending  on  the  boundary  condition) 
to  separate  a  and  b  .  Finally  the  frequencies  required 

would  be  quite  low,  thus  perhaps  complicating  the  experimental 
task  of  frequency  analysis.  Further  study  is  being  given  to 
this  problem. 
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